Spin-glass state of vortices in YBa2Cu30y and La2- 
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Highly disordered magnetism confined to individual weakly interacting vortices is detected by 
muon spin rotation (/xSR) in two different families of high-transition-temperature (Tc) supercon- 
ductors, but only in samples on the low-doping side of the low-temperature normal state metal-to- 
insulator crossover (MIC). The results support an extended quantum phase transition (QPT) theory 
of competing magnetic and superconducting orders that incorporates the coupling between Cu02 
planes. Contrary to what has been inferred from previous experiments, the static magnetism that 
coexists with superconductivity near the field-induced QPT is not ordered. Our findings unravel the 
mystery of the MIC and establish that the normal state of high- Tc superconductors is ubiquitously 
governed by a magnetic quantum critical point (QCP) in the superconducting phase. 

PACS numbers: 74.72.-h, 74.25.Ha, 74.25.Qt, 76.75. 



I. INTRODUCTION 

For two decades, arrival at a firm theory for high-Tc 
superconductivity has been hindered by an incomplete 
characterization of the phase diagram for cupratc mate- 
rials. In zero field, /iSR^i^i^iii^i^ NMR/NQRl and neu- 
tron scattering^ studies show that static (or quasistatic) 
magnetism coexists with superconductivity in the under- 
doped regime. Field-induced or enhanced static mag- 
netic order has also been clearly detected in underdoped 
La2_^Sr^Cu04 {hSCO:x) ^^^^ Pri_^LaCe^Cu04?i^>i^ 
and La2Cu04+j/^ii^ by neutron scattering, and in under- 
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doped Pr2-xCe:,Cu04 (PCCO),!^ Pri_:rLaCexCu04 
and LSCOiceiS, by /iSR. The neutron studies on LSCOia; 
and La2Cu04+j^ support a proposed phase diagram by 
Demler et al^ in which the pure superconductor under- 
goes a QPT to a phase of coexisting static magnetic and 
superconducting orders. A similar phase transition com- 
patible with the theory of Ref.^5. has also been observed 
in CeRhlns.^^ Still the general applicability of this QPT 
model is questionable, since field-induced static magnetic 
order has not been established in any of the other hole- 
doped cuprates. 

An important detail in the model of Refi^ is the as- 
sumption that the vortices are two-dimensional (2D). In 
this case the competing order is stabilized only when 
there is strong overlap of the 2D vortices within a Cu02 
layer. When this happens long-range magnetic order is 
established. However, Lake et al^ have shown that 
the field-induced order in LSCOiO.lO is in fact three- 
dimensional (3D), implying significant interlayer cou- 
pling. Furthermore, a neutron//^SR study of LSCOiO.lO 
concluded that the vortices themselves are SD.^'^ Follow- 
ing the work of Ref,-^, Kivelson et al. showed that com- 
peting order can be stabilized about a nearly isolated 
3D vortex,^ The field-induced QPT in this extended 3D 



model is argued to be to a coexistence phase in which 
the spatial dependence of the competing order is sub- 
stantially non-uniform. 

Here we show that there is a generic field-induced 
transition to a coexistence phase where spin-glass-like 
(SG) magnetism is confined to weakly interacting 3D 
vortices. The detection of this phase implies that the 
QPT previously identified in LSCOia; by neutrons scat- 
tering is simply a crossover to a situation where com- 
peting static magnetism is spatially uniform in the sam- 
ple. Furthermore, we identify the "true" field-induced 
QPT as occurring near the critical doping for the low- 
temperature normal-state MIC that occurs at a non- 
universal doping concentration in cuprate superconduc- 
tors. The insulating side of the normal-state MIC is 
characterized by a log(l/T) divergence of the in-plane 
resistivity p„f,^i2^i^i2^i2ai^i^i^ but it has also been 
indirectly identified by electronic thermal conductivity 
measurementsj^^i^i^ 



II. EXPERIMENT 

Muon spin rotation/relaxation (/xSR) measurements 
were performed at TRIUMF, Canada on LSCOra: and 
YBa2Cu30y (YBCO:y) single crystals on either side 
of the previously determined critical dopings Xc ~ 
0.16 (Refs^6'33.34) and ~ 6.55 (Ref.^S) for the low- 
temperature MIC. The LSCOix single crystals were 
grown by the traveling-solvent floating-zone technique;^ 
whereas the YBCOiy single crystals were grown by a self- 
flux method in fabricated BaZrOs cruciblesi^ 

The /iSR method involves the implantation of nearly 
100 % spin polarized positive muons into the sample. 
Like a tiny bar magnet, the magnetic moment of the 
muon precesses about the local magnetic field B with an 
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FIG. 1: (Color online) ZF-/iSR time spectra at different 
temperatures for (a) YBCO:6.67, (b) YBCO:6.57 and (c) 
YBCO:6.50. In all cases the initial muon spin polarization 
P(0) was parallel to the c-axis of the single crystals. The 
solid curves through the data points are fits described in the 
main text. 



angular frequency =7^i?, where 7^ = 0.0852 /is^^ G^^ 
is the muon gyromagnetic ratio. By measuring the 
time evolution of the polarization of the ensemble of 
muon spins P{t) via the anisotropic distribution of decay 
positrons, the internal magnetic field distribution n{B) 
of the sample is determined.- - As described below, /iSR 
measurements were first carried out in zero external field 
to search for static electronic moments. The vortex cores 
were then probed by transverse-field (TF) /iSR, with the 
applied magnetic field perpendicular to the Cu02 layers. 



III. ZERO-FIELD MEASUREMENTS 

Figures [T] and m show ZF-/xSR time spectra for some of 
the samples. Defining the direction of the initial muon 
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FIG. 2: (Color online) ZF-/iSR time spectra at different tem- 
peratures for the lowest doped samples (a) YBCO:6.46, and 
(b) LSCO:0.145, recorded with the initial muon spin polar- 
ization P(0) perpendicular to the c-axis of the single crystals. 
The solid curves through the data points are fits described in 
the main text. 

spin polarization P(0) to be parallel to the z-axis, the 
ZF-fiSK or "asymmetry" spectrum has the form 

A{t)=aoP,{t)=aoG,{t), (1) 

where ag is the initial asymmetry and Gz{t) is a relax- 
ation function, fn all cases the spectra are well described 
by the following ZF relaxation function 

G,(i) = GrWexp(-AO, (2) 

where 

GrW = J + ^(l-A^i^)exp(-iA^t^) , (3) 

is the static Gaussian Kubo-Toyabe function (KT) typ- 
ically used to describe relaxation due to nuclear dipole 
fields, and A is an additional exponential relaxation rate, 
fn the absence of static or slowly fluctuating electronic 
moments, the relaxation of the ZF-/iSR signal is caused 
solely by the nuclear dipoles. fn this case the relax- 
ation is expected to be independent of temperature, as 
observed in Figs. [1] and [2l Fitted values for A and A 
are given in Table [D The measurements on YBCO:6.46 
and LSCO:0.145 were done using a different spectrom- 
eter, and with the initial muon spin polarization P(0) 
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FIG. 3: (Color online) Doping dependence of the ^SR line shapes for YBa2Cu30y at H = 5 kOe and T — 2.5 K. (a) Full fiSH line 
shape for samples above {y — 6.57) and below (y — GAG) the critical doping j/c = 6.55 for the MIC. (b), (c), (d), (e) Blowups of 
the 'tail' region of the fiSH line shapes above (j/ = 6.67, 6.60 and 6.57) and below {y — 6.50 and 6.46) yc = 6.55. For comparison, 
all line shapes have been normalized as described in the main text. 



TABLE I: Resuhs of fits of the ZF-/iSR time spectra to 
Eq. ((2)| . The fits are shown as solid curves in Figs. [1] and 




Sample A (/is "'") A (/is ^) Polarization 

YBCO:6.50 0.1120(3) 0.0381 P(0) || c 

YBCO:6.57 0.1194(2) 0.0244 P(0) || c 

YBCO:6.67 0.1207(1) 0.0332 P(0) || c 

YBCO:6.46 0.1277 0.1044 P(0) _L c 

LSCO:0.145 0.195(10) 0.11(2) P(0) ±c 



perpendicular, rather than parallel to the c-axis. In 
this geometry the relaxation rate is larger due to the 
anisotropy of the muon-nuclear dipole interaction.'^^ The 
hole doping dependence of A in YBCO:?/ is explained 
by a change in the ratio of muons stopping near the 
0(1) and 0(4) oxygen sites4S While there is a resid- 
ual exponential relaxation rate for all samples, A is in- 
dependent of both temperature and hole doping con- 
centration. Thus there is no evidence from the ZF- 
/iSR spectra for static electronic moments in our sam- 
ples, which is an essential requirement for establishing 
the presence of hidden competing magnetic order. We 
remark that the temperature-independent exponential 
component may come from the fraction of muons miss- 
ing the sample and avoiding the background suppression 
scheme of the spectrometer. Furthermore, the measure- 
ments here do not rule out the presence of a weak tem- 
perature dependent relaxation rate found in earlier high 
precision ZF-/iSR measurements of YBCO:y42ii^ 



IV. TRANSVERSE-FIELD MEASUREMENTS 

In a transverse field, the muon spin precesses in a plane 
perpendicular to the field axis. In this case the asymme- 
try spectrum is 

A{t) = aoP^ (t) = aoG, {t) cos{-f^Bt) , (4) 

where Gx (t) is the transverse muon spin relaxation func- 
tion and B is the local field at the muon site. In the vor- 
tex state the internal magnetic field is spatially inhomo- 
geneous, and the TF-/xSR signal for a perfectly ordered 
flux-line lattice (FLL) is described by the polarization 
function 

P,(t) =^cos[7^B(r,)i], (5) 

i 

where the sum is over all sites in the real-space unit cell 
of the FLL and B{ri) is the local field at position = 
{xi,yi). A Fourier transform of Px(t) 

/•OC 

n{B)^ P^(t)e-*(''^-^*'e-'^AtV2^5^ (g) 
Jo 

often called the '/iSR line shape', provides a fairly ac- 
curate visual illustration of the internal magnetic field 
distribution sensed by the muons. Here exp{-alt^ /2) 
is a Gaussian apodization function used to suppress the 
"ringing" effect of the finite time range of Px{t). Fig- 
ure El^a) shows a couple of examples of the /iSR line 
shape for YBCO:y at H ^ 5 kOe and T = 2.5 K. The 
asymmetric line shape for YBCO:6.57 is typical of the 
field distribution for a 3D FLL.^® Specifically, the 'high- 
field' tail corresponds to the spatial region in and around 
the vortex cores. 
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FIG. 4: (Color online) Doping, temperature and magnetic field dependences of the ^SK line shapes for La2-iSra:Cu04. The 
center panel shows the full /xSR line shapes at H — 5 kOe and T = 5 K, above (j: = 0.176 and a; = 0.166) and below (a:: = 0.145) 
the critical doping a::c = 0.16 for the MIC. Panels (a)-(c) show the field dependence of the 'tail' regions of the normalized /xSR 
line shapes. Panels (d)-(f) show the temperature dependence of the 'tail' regions of the normalized /xSR line shapes. 



Figures [3fb)-(e) show blowups of the 'tail' regions of 
the Fourier transforms of TF-/xSR spectra measured in 
the vortex state of YBCO.y near the critical doping 
He ~ 6.55. For comparison, the line shapes have been 
normalized to their respective peak amplitude ripk(-B). 
Furthermore, to account for differences in the in-plane 
magnetic penetration depth, the widths of the line shapes 
have been made equivalent by rescaling the horizontal 
B—Bq axis, where Bq is the applied magnetic field. Above 
Uc the fiSK line shapes for ?/ = 6.67, j/ = 6.60 and y = 6.57 
are identical. However, at y = 6.50 there is a clear sup- 
pression of the high-field tail, corresponding to the spa- 
tial region of the vortex cores. Note that in a previous 
high-field study of YBCO:6.50 the high-field tail was ar- 
gued to be enhanced rather than suppressed42. However, 
this conclusion was based on a comparison of the /iSR 
line shape to an assumed theoretical curve for n{B). As 
observed in Fig. Ele), at ?/ = 6.46 the suppression of the 
high-field tail is accompanied by the appearance of a low- 
field tail. 

As shown in Fig. 31 similar differences are observed 
between the /iSR line shapes of LSCOia; above and below 
the critical doping a;c = 0.16 for the MIC. With increasing 
magnetic field the differences between the tails of the 
line shapes are enhanced [see Figs Sl^a), |3Jb) and (He)]. 



On the other hand, with increasing temperature the /iSR 
line shape of LSCO:0.145 becomes more like that of the 
samples above Xc — 0.16 [see Figs. IHJd), IDJe) and|D^f)]. 
In the next section we explain how field-induced static 
electronic moments in the samples on the low-doping side 
of the MIC accounts for both the suppression of the high- 
field tail and the appearance of a low-field tail. 

V. DATA ANALYSIS 

While the change in the local magnetic field distribu- 
tion in the region of the vortex cores is evident from a 
visual inspection of the /iSR line shapes in Figs. [3] and 
21 it is constructive to consider a simple analysis of the 
TF-/iSR time spectra. Recently we carried out a compre- 
hensive analysis of the /xSR line shape in y > 6.57 single 
crystals4^ There we showed that the TF-/iSR signal is 
well described by the polarization function 

P.{t) = e--'«*'/2 cos[7^i3(r0t] , (7) 

i 

where the Gaussian function exp(— accounts for 
additional relaxation due to FLL disorder and nuclear 
dipole moments, the sum is over all sites in an hexagonal 
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FLL, and B{ri) is the following analytical solution of the 
Ginzburg-Landau equations'*'^ 

Here G are the reciprocal lattice vectors of the FLL, Bq is 
the average internal magnetic field, F{G) = uKi{u) is a 
cutoff function for the G sum, Ki [u] is a modified Bessel 
function, and u — \/2£_abG. The cutoff function F{G) de- 
pends on the spatial profile of the superconducting order 
parameter at the center of the vortex core. Consequently, 
the parameter ^ab is a measure of the vortex core size. As 
explained in Ref.— , only the H^O extrapolated value of 
Xab is a true measure of the magnetic penetration depth, 
since at finite H this parameter absorbs deviations of 
B{ri) from Eq. dH). 

In Ref.-^ it was assumed that the unusual /iSR line 
shape of YBCO:6.50 results from static antiferromag- 
netic order in the vortex cores. However, field-induced 
static magnetic order has never been observed in YBCO:y 
by neutron scattering. Furthermore, Khaykovich et al. 
have shown by neutron scattering that static magnetic 
order occurs in LSCO:0.144 only above 7J « 30 kOeJ^ 
Thus the AiSR line shapes of YBCO:6.46, YBCO:6.50 
and LSCO:0.145 presented here for if < 15 kOe cannot 
be explained by static magnetic order in and around the 
vortex cores. Instead we consider the possibility that 
the weak fields considered here induce disordered static 
magnetism, which is not ruled out by the neutron scatter- 
ing experiments. A polarization function that describes 
the case of disordered static electronic moments in and 
around the vortex cores is 

P.it) = e-^'«*'/2^exp(-Ae-('-'/«-)'0cos[7^S(r,)t]. 

(9) 

This equation simply says that a muon stopping at posi- 
tion Ti in the FLL experiences a Lorentzian distribution 
of fields typical of a SG system that results in an ex- 
ponential decay of P{t). Furthermore, the exponential 
relaxation rate A, and hence the width of the field dis- 
tribution, is assumed to decrease with increased distance 
from the center of the vortex core. 

Figure [5] shows the real part of the Fourier transform 

/■oc 

n{B)= P^(i)e-''^^-«*e-'"A*'/^dt, (10) 

where [t) is calculated from Eq. ^ for the case cFcS = 0. 
The vortex cores are non-magnetic for the case A = 0. 
When A is non-zero the high-field tail of the line shape 
is suppressed. With increasing A the high-field tail is 
further suppressed, and a low-field tail develops. While 
the change in the high-field tail is most recognizable, the 
appearance of the low-field tail depends on the width of 
the SG Lorentzian field distribution relative to the line 
width of n{B) for the FLL. For example, in Figs. EKb) 
and the Lorentzian field distribution of the static 




B - B„ (G) 

FIG. 5; (Color online) Real part of the Fourier transform of 
Eq. (|9]) calculated for two values of A (0 /is~^ and 5 Ats^^), 
with (Teff = and (a) A„6 = 2000 A, = 50 A, and Bo = 
5 kG, (b) Aab = 2000 A, Cab = 100 A, and Bo = 5 kG, and (c) 
Aab = 3000 A, Cab = 50 A, and Bo = 15 kG. All of the Fourier 
transforms have been generated with a Gaussian apodization 
of width (TA = 0.2 /is~^. The insets show blowups of the 
bottom portion of the same Fourier transforms. 



magnetism is broad enough to extend beyond the low- 
field cutoff of the field distribution of the FLL. 

Figure [6] shows the temperature dependence of A ob- 
tained from fits of the TF-/iSR signal to Eq. ^ for 
some of the samples. Since A « for LSCO:0.176 and 
YBCO:6.60, we conclude that the vortex cores are free 
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FIG. 6: (Color online) Temperature dependence of A from 
fits of the TF-/iSR time spectra at H = 5 kOe to Eq. ^ for 
LSCOia; and YBCOij/ samples above and below the critical 
dopings Xc = 0.16 and j/c ~ 6.55, respectively. At the lowest 
temperature the fits yield \ab ~ 2505(30) A, £,ab ~ 34(3) A, 
and A = 6.5(1.9) /is'^ for LSCO:0.145, and A,,^ = 1950(56) A, 
Cab = 90(10) A, and A = 5.9(1.8) fis'^ for YBCO:6.50. Note 
that the vortex core size in YBCO:y at low field is large due 
to proximity-induced superconductivity on the CuO chain 
layers.— 



of static magnetism. On the other hand, the diverg- 
ing temperature dependence of A for LSCO:0.145 and 
YBCO:6.50 indicates a static broadening of the internal 
magnetic field distribution associated with the spatial re- 
gion of the vortex cores. In other words, the increase in A 
with decreasing temperature is consistent with a slowing 
down of fluctuating Cu spins. Moreover, since the value 
of A does not saturate down to T — 2.5 K, the tempera- 
ture dependence of A is consistent with an approach to 
a second-order magnetic phase transition at T = K. 

The main changes in the /iSR line shape across the 
critical doping for the low-temperature MIC are now 
understandable. From the fitted values of A, the half- 
width at half-maximum of the Lorentzian field distri- 
bution assumed in Eq. ^ is approximately ±70 G in 
both LSCO:0.145 and YBCO:6.50 at T = 2.5 K. For 
LSCO:0.145 this is broad enough to affect the low-field 
tail. With increasing field, the density of magnetic vor- 
tices increases, while the field inhomogeneity of the FLL 
decreases. Consequently, at higher magnetic field the 
static broadening of the /iSR line shape by the magnetism 
becomes more discernable [see Figs. |ll[a)-(c)]. With in- 
creasing temperature, the simultaneous loss of the low- 
field tail and the recovery of the high-field tail of the 
LSCO:0145 fine shape [see Figs. Ill^d)-(f)] signifies ther- 
mal destruction of the static magnetism in and around 
the vortex cores. 
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FIG. 7: (Color online) Low-temperature ZF-^SRtime spectra 
for YBCO:y on either side of the critical doping for the MIC 
at j/c = 6.55. 



VI. DISCUSSION 

A. Dimensional crossover? 

While the onset of SG magnetism in and around the 
vortex cores fully explains our experimental observations, 
we note that the change in the /xSR line shape across 
the MIC is somewhat reminiscent of that observed across 
the 3D-to-2D vortex crossover field in highly anisotropic 
Bi2+a;Sr2_a;CaCu208+5 (BSCC0)4^ In this case random 
pinning-induced misalignment of the stacked 2D 'pan- 
cake' vortices that comprise the 3D fiux lines in BSSCO, 
narrows and reduces the asymmetry of the /iSR line 
shape. The two ingredients necessary for such a crossover 
are weak coupling between the Cu02 planes and a source 
of pinning. 

If the vortices in LSCO:0.145 are quasi-2D, the grad- 
ual recovery of an asymmetric line shape at higher T 
that is observed in Fig. [4] signifies thermal depinning 
of the vortices and a return to an ordered 3D FLL. 
Such a scenario has been observed in BSCC04^ How- 
ever, as already mentioned, the vortices in LSCO:0.10 are 
known to be 3D.— Since the effective mass anisotropy 
T = ■\/m*/m*^ increases with decreasing hole doping 
concentration ; a novel mechanism that softens the 
vortex lines at higher doping would be needed to explain 
the LSCO:G.145 line shapes. As for YBCO:y, mutual 
inductance measurements show that even severely un- 
derdoped samples are quasi-2D only near Tc4^ The weak 
field dependence of the Josephson plasma resonance in 
YBCO:6.50 at low T is also consistent with 3D vortices. — 
Thus the extreme vortex anisotropy necessary for a 2D- 
to-3D crossover does not seem to occur at the hole-doping 
concentrations investigated here. 

Assuming this were not the case, one could imag- 
ine an abrupt onset of disorder at the MIC that drives 
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FIG. 8: (Color online) Schematic zero-temperature H versus hole doping concentration p phase diagram for (a) 2D vortices 
(adapted from figure 1 of Ref.— ), and (b) 3D vortices (adapted from figure 1 of Ref.— ). ^ and <^ denote the expectation values 
of the superconducting and competing order parameters, respectively. In each panel the solid blue curve is a QPT and the solid 
blue dot denoted Pc is a QCP. Note that the QPT in (a) becomes a crossover (red dashed curve) in (b). Below the crossover, 
the competing order is spatially non-uniform. In addition, the — > extrapolation of the QPT in (b) at p\ is an 'avoided' 
QCP (open blue circle). Note the Meissner phase is not shown in (a) or (b). 



both the logarithmic divergence of the normal-state re- 
sistivity and a 3D-to-2D vortex crossover. Since the lo- 
cation of the MIC in LSCO:x has been independently 
confirmed,— i^i^ such pinning would have to be intrinsic 
to the material. However, disorder in LSCO:a; decreases 
with decreasing Sr doping, and disorder due to excess 
or deficient oxygen primarily affects the lightly-doped 
and overdoped regimes, respectively. Likewise, the onset 
of pinning below the MIC is inconsistent with ortho-II 
ordering in YBCO:6.50, which reduces random pinning 
by oxygen disorder and defects. An abrupt redistribu- 
tion of charge at the MIC is also not supported by our 
own ZF-/iSR measurements. Previously, ZF-/iSR studies 
of cuprate superconductors have demonstrated a sensi- 
tivity to charge-poor magnetic regions^, and to charge 
correlationsiS. However, as shown in Fig. [3 the low- 
temperature ZF-/uSR signal does not change across the 
MIC. 



B. Avoided quantum criticality 

Figure [8] shows a comparison of the proposed phase di- 
agrams for competing order in the cuprates for the case 
of 2D (Ref.20) and 3D (Ref.^^) vortices. The major dif- 
ference is that the inclusion of the interlayer coupling 
allows the competing phase to be stabilized in nearly iso- 
lated vortices, thus altering the position and character 
of the QPT. There are two key predictions of the ex- 
tended theory for 3D vortices that are confirmed by our 
experiments. The first is that there exists a coexistence 
phase of spatially inhomogeneous competing order. At 
low fields where the interaction between vortex lines is 
weaker, we have detected SG magnetism that is local- 



ized in and around the vortex cores. This means that 
the competing order initially stabilized at the QPT is 
not static magnetic order as previously established, but 
rather is highly disordered static magnetism. The com- 
peting magnetism is characterized by a local order pa- 
rameter, namely, the mean squared local magnetization. 
With increasing field, stronger overlap of the magnetism 
around neighboring vortices may lead to a co-operative 
bulk crossover to long-range magnetic order, as is ap- 
parently the case in Lai.856Sro.i44Cu04.— While field- 
induced static magnetic order has not been detected in 
YBCO:y, it is worth noting that the MIC occurs at a 
much lower hole doping concentration (pc ~ 0.1) than 
in LSCO:a; (pc = 0.16). Consequently, very high mag- 
netic fields are likely needed to induce long-range mag- 
netic order in YBCO:?/. In contrast, very weak fields 
were shown to induce magnetic order in PCCO samples 
that are below the MIC crossover at pc~0.16iii This is 
understandable, since in zero field the superconducting 
phase of PCCO is in close proximity to the pure antifer- 
romagnetic phase where the competing order parameter 
is spatially uniform throughout the sample. 

Another key prediction of Refi^ is that there is an 
'avoided' QCP at H — 0, meaning that the QCP lies at 
a lower doping than the extrapolated H^Q value of the 
field-induced QPT. This is shown in Fig. [Sfb). Consis- 
tent with this idea, ZF-/zSR studies of pure LSC0:a4ii 
and YBCOyijS- indicate that the onset temperature for 
coexisting static magnetism and superconductivity ex- 
trapolates to zero below the critical doping for the MIC. 
While this is well below the doping concentration p = 0.19 
that Tallon and others^i^ have advocated to be a univer- 
sal QCP in the cuprates, we stress that our study does 
not prohibit the existence of more than one QCP under 
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FIG. 9: (Color online) Schematic r = K phase diagram de- 
duced from this study. The normal and superconducting (SC) 
phases occur above and below the upper critical field Hc2{p), 
respectively. The solid vertical curve at Pc{H) is a QPT coin- 
ciding with the low-T normal-state MIC. Below Hc2{p), Pc{H) 
separates a pure SC phase from a SC phase with coexisting 
static magnetism (SM). Immediately to the left of Pc{H) the 
SM is disordered, becoming spatially uniform (and possibly 
ordered) above the dashed curve. The open circle is the pre- 
dicted 'avoided' QCP,™ whereas the solid circle indicates the 
'true' QCP at H = 0. 

the superconducting 'dome'. 

VII. CONCLUSIONS 

Our experiments clearly demonstrate a change in the 
internal magnetic field distribution of the vortex state 



across the critical doping concentration for the low- 
temperature MIC in two hole-doped high- Tc supercon- 
ductors. We have shown that the occurrence of SG mag- 
netism in and around weakly interacting vortices is the 
most likely source of the observed changes. In Fig. [9] 
we show a generic zero-temperature phase diagram that 
is compatible with the present and previous works. We 
conclude that the strange localization of charge below 
the MIC stems from competing static magnetism that is 
stabilized when superconductivity is suppressed by the 
applied field. While others have hypothesized that mag- 
netism is the cause of the peculiar localization of charge, 
the experiments here establish that static magnetism not 
present in zero external field does appear in an applied 
magnetic field immediately below the critical doping for 
the MIC. Magnetism plays a prominent role in at least 
one theory for the MIC. In particular, Marchetti et al. 
have used a spin-charge gauge approach to show that the 
MIC can arise from a competition between short-range 
magnetic order and the dissipative motion of the charge 
carriers/S The experiments here do not rule out the pos- 
sibility that there are short-range spin correlations in the 
field-induced magnetism immediately below the MIC. 
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